Probing nonperturbative QED with optimally focused laser pulses 
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We study nonperturbative pair production in intense, focused laser fields called 'e-dipole' pulses. 
We address the conditions required, such as the quality of the vacuum, for reaching high intensities 
without initiating beam-depleting cascades, the number of pairs which can be created, and exper- 
imental detection of the created pairs. We find that e-dipole pulses ofi^er an optimal method of 
investigating nonperturbative QED. 



Understanding quantum field theory in the nonpertur- 
bative regime remains a challenging theoretical and ex- 
perimental issue. Recent advances in technology have 
spurred interest in the possibility of using intense lasers 
to probe quantum vacuum phenomena, including nonper- 
turbative electron-positron pair production [U [2] . Such 
studies are complementary to those of, e.g., heavy ion 
collisions, with the possible advantage that laser systems 
provide a more controlled environment. 

It is well known that nonperturbative pair creation 
is strongly suppressed below the Sauter-Schwinger limit 
Es = m^c'leh ~ IQi* V/m [3 g], and while field- 
strengths of this scale are typical of QED, creating them 
on a macroscopic (laboratory) scale will remain out of 
our reach for the foreseeable future. Various mecha- 
nisms for stimulating pair production have therefore been 
proposed; these include colliding an intense laser pulse 
with high energy photons [5] , electrons [B] , or other laser 
pulses . The number of pairs which can be created 
in the collision of laser pulses is sensitive not just to the 
field amplitude [10], but also to the field structure [11]. 
The same applies to particle dynamics in general. In 
simulating intense laser-matter interactions it is there- 
fore necessary to employ realistic field models. 

It was pointed out in [H], see also jlSj, that there 
is a potential obstacle to reaching the Sauter-Schwinger 
limit, or high intensities in general. In an experiment, 
the presence of stray particles resulting from imperfect 
vacuum, or the generation of hard radiation, could trig- 
ger an avalanche of particle production. This pertur- 
batively triggered cascade could hinder any observation 
of nonperturbative effects, both via beam depletion and 
by producing a background which swamps the signals of 
interest. Hence, it is important to understand the con- 
ditions leading to cascades, and what backgrounds they 
produce. 

In this paper we consider Schwinger pair production 
in 'e-dipole' pulses [14]. These are exact, closed form 
solutions of Maxwell's equations in vacuum. They exhibit 
optimal focusing efficiency, i.e. they produce the highest 



possible peak field-strength for a given input power [T5] 
and describe genuinely pulsed fields, having finite energy 
and finite extent in all four spacetime directions. We will 
show in a separate article that it is possible, using a new 
focusing system, to create a close approximation to an e- 
dipole pulse, hence pushing laser systems to near-optimal 
conditions [T^j. Here we focus on exact e-dipole pulses 
in order to provide analytic estimates; we will show that 
the structure of e-dipole pulses may reduce the chances 
of cascade formation, thus allowing laser experiments to 
approach higher intensity regimes. 

This paper is organised as follows. We begin by in- 
troducing e-dipole pulses. We then analyse particle mo- 
tion in such pulses, showing that stray particles in an 
experiment are not driven into the high intensity part 
of the field, provided one begins with a vacuum of suffi- 
cient quality. We estimate the level of vacuum required 
to both keep particles away from the focus, and to limit 
the generation of hard photons into the focus. We then 
calculate the number of pairs which could be produced 
by our pulse and compare with other pulse models in 
the literature. The behaviour of electron-positron pairs 
post-creation is then analysed before we conclude. 

Dipole pulses:- There is a limit to the focusing effi- 
ciency of a given laser system. Consider first an unfo- 
cused, broad laser pulse which can be regarded as nearly 
monochromatic. The peak field strength Eq which can be 
obtained by focusing monochromatic light of wavelength 
A and cycle-averaged power P is bounded above; defining 
P = P/(47reoc), the peak field obeys [15] 
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The focused fields which saturate this upper bound are 
called e-dipole pulses, see [2] for full details. i?-dipole 
pulses take their name from structural similarities with 
dipole fields, but do not contain singularities. They de- 
scribe a converging pulse of light, with the ideal case of 
47r focusing. (The formation of this focused pulse can be 
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pictured as the reverse process of emission from a dipole 
antenna.) The e-dipole field has the following form. Let 
— + ^ and define the vector Z by 

Z^z^[g{t + R/c)~g{t~R/c)\ , (2) 

in which the 'driving function' g is arbitrary and the 'vir- 
tual dipole moment' d is a constant. (Using the fact 
that e-dipole pulse can be generated by focusing laser 
fields [13], we will shortly relate g and d to input laser 
parameters.) The electromagnetic fields of an e-dipole 
pulse are then given in terms of Z by 



E = -VxVxZ, B=:-4vxZ 



(3) 



These fields are exact solutions to Maxwell's equations in 
vacuum. It is easily confirmed that there is no singularity 
at i? = 0, the focus point. Far from the focus, the electric 
and magnetic field amplitudes of the e-dipole pulse are 
proportional to dR~^g(t ± and have an angular dis- 
tribution proportional to sin^ 9, with 9 the polar angle. 
In the focus one has E(0,t) = z^'g{t), B(0,0 = 0. In 
general we are interested in pulses with, say, Gaussian 
frequency spreads as for the following driving function 



g(T) = e-(^'/^')'"4sin(c^r) 



(4) 



in which lj is the central frequency and D is the full- 
width-half-maximum duration (i.e. the intensity ~ 
drops to 1/2 its peak value at i = ±D/2). We refer to 
such a pulse as 'quasi-Gaussian' since, far from the focus, 
the envelope is g, which has the same Gaussian frequency 
spread as g. It follows from ([3| that such pulses are 
compactly supported in all four space-time dimensions, 
which is one of the advantages of e-dipole pulses over 
other field models in the literature; they describe genuine 
pulses, without sacrificing Maxwell's equations. 

The virtual dipole moment d can be expressed as 
a function of input energy, or power, by using energy 
conservation in the far-field. To illustrate, consider a 
monochromatic driving field, wavelength A and cycle- 
averaged power P := ^Ppeafc- The peak focused field 
strength _Bo, incoming power, and d are then related by 
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Note that ([s]) is optimal, saturating the upper bound in 
([T]). Similar relations can be evaluated for the quasi- 
Gaussian pulse [14 . This leads to a simple measure for 
comparing different pulse models. We define the focus- 
ing efficiency parameter Je via Eq = Je^] this relates 
the incoming power P to peak field-strength Eq. In or- 
der to compare with other models, this definition applies 
to fields which are monochromatic in the far field re- 
gion. In Table [l] we compare the focusing parameters for 
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TABLE I. Focussing efficiency /_e (input power to peak elec- 
tric field) for three beam configurations. We use A = 0.1 to 
compare with the literature (A ^ 1 is required jTU])- The 
e-dipole pulse is most efficient. 



three fields: the theoretical maximum from Fedotov's so- 
lution |T7], two colliding Narozhny-Fofanov beams [TS] 
and the e-dipole field (with monochromatic g). The e- 
dipole pulse is most efficient. 

Given a particular pulse model, one can increase (de- 
crease) the focal electric (magnetic) field-strength by col- 
liding several such pulses [S]. This increases the total 
focusing efficiency and is why, in Table [l| we compare 
with two colliding beams from |18j , as advocated for pair 
production in [19]. Using two pulses eliminates the focal 
magnetic field of the single pulse solution, and increases 
the peak field strength by a factor of The poten- 
tial increase using this method is limited, though; the 
behaviour of the efficiency depends strongly on the fo- 
cusing geometry, due to the requirement that the total 
angle of incoming radiation cannot exceed 47r. As multi- 
ple beams start to overlap, they interfere, and the peak 
field strength is reduced. In particular, the 'y^ increase 
in field strength from n pulses' is limited, for fixed and 
reasonable focusing angles, to rather small n. 

Interaction with particles within the chamber:- A po- 
tential barrier to the observation of nonperturbative 
pair production is, simply, perturbative pair production, 
which is well understood. This and other processes can 
be initiated by particles within the imperfect vacuum of a 
chamber prepared for a laser experiment. Due to the ef- 
fect of the ponderomotive force, an incoming laser pulse 
can 'grab' stray particles and accelerate them into the 
focus, leading to the emission of high energy photons 
which create pairs, which emit photons which create fur- 
ther pairs. The products of this cascade can obscure any 
nonperturbatively produced pairs, or the resulting beam 
depletion can be so great as to prevent sufficiently high 
intensities for Schwinger pair creation being reached jl2| . 
If, though, we could guarantee that within the focal vol- 
ume there were no high energy electrons and photons ini- 
tiating perturbative processes (or rather, that such events 
had a low probability), then we could focus beams to 
the intensities required to trigger nonperturbative pair 
creation, without first initiating cascades, which is the 
experimental result of interest. 

To investigate this prospect, one can determine the vol- 
ume of space from which electrons can initiate cascades, 
as a function of the initial electron density (initial level 
of vacuum), and then require that this volume contain 
less than one particle. We therefore carried out a simula- 
tion of a large number of initially uniformly distributed 
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FIG. 1. Simulation results for electron motion under the 
influence of an e-dipole pulse. Left scale: initial number (No) 
and final average number (A^) of electrons dragged into a 
sphere of radius R about the focus, as a function of R. Ini- 
tial particle density 10^ cm"''. Right scale: maximum values 
of X fi'iid ^, taken over all electrons outside the sphere with 
radius R. 



electrons moving (with radiation reaction accounted for 
by the Landau-Lifshitz equation [20]) in the fields of 
a converging e-dipole pulse, wavelength 810 nm, pulse 
duration 30 fs (consistent with candidate Ti-sapphire 
technology for future laser facilities [211 [22) and total 
power P = 1000 PW (as is estimated, below, to be re- 
quired for pair creation, see also Table |ll| . Laser experi- 
ments are now routinely performed in technical vacua of 
stray molecule density 10~®-10~^mbar, and lower pres- 
sures are accessible. We will assume an initial density of 
lO^cm"'^, equivalent to a pressure of lO^^^mbar. This 
figure is at the limit of what is achievable today, so it a 
reasonable assumption for future facilities. The incoming 
laser pulse will ionise stray molecules and produce elec- 
trons. Now, we are interested only in whether electrons 
are pushed into the focus or not, based on their initial 
position. Thus, for the determination of the required 
vacuum quality, we can associate every molecule with a 
single electron, independent of the level of ionisation. We 
therefore take the initial electron density to be lO^cm"'^. 

The results of our simulation are shown in Fig. [ij We 
show here the average number of electrons initially dis- 
tributed within (A^Oi calculated from the initial density), 
or dragged into (N), a sphere of radius R around the fo- 
cus, as a function of R. We see that for the chosen initial 
density, the final number N of particles drops below one 
at a radius of around R* = BO/im; this is much larger 
than the focus of the dipolc pulse. 

We immediately perform two checks on this result. 
First, we should only base our conclusions on data from 
regions of space in which quantum effects are negligible. 
Second, though particles remain outside the focus, they 
may still be accelerated such that they emit hard pho- 



tons (capable of initiating cascades). In order to verify 
our result, we tracked in our simulation two parameters, 
X and The first, Xi is the quantum efficiency param- 
eter of the particles, calculated along their orbits [23j . 
This parameter estimates the importance of quantum ef- 
fects and is, for a particle of momentum in a field with 
energy-momentum tensor T^^, |24j 
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In the electron's instantaneous rest frame, x is the ratio 
of the energy seen by the electron, to the electron's rest 
energy. As x approaches unity, quantum effects start to 
become important; in high intensity fields x is the in- 
put parameter for determining the probability of QED 
processes [53] ■ The second parameter, ^, estimates the 
energy of emitted photons, in ratio to twice the rest mass 
of an electron. Recall that the typical energy hiOs of pho- 
tons emitted by an ultra-relativistic electron is given in 
the synchrotron approximation by Ws = 7^(e_E/mc) [20j . 
We therefore define 
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So, for ^ ^ 1 we should be able to neglect the emission 
of hard photons with energy approaching 2mc^ . 

Returning to our results. Fig. [T] also shows the max- 
imum value (over all electrons) of x and ^ outside the 
sphere of radius R. We can conclude first that the clas- 
sical description is valid (x ^ 1) at least outside the 
sphere of radius R* . Hence, our estimate of the density 
increase within this sphere (due to dragging from outside 
it) can be trusted. Second, since ^ < 1, particles outside 
the sphere do not generate photons suitable for initiating 
pair creation within the sphere. Hence we can conclude 
from Fig. [ijthat an initial density of less than 10^ cm~^ 
leads to an average of well less than one particle existing 
in or entering into the sphere of radius R* , and an ab- 
sence of photons sufficiently hard to create pairs within 
the sphere. 

Our estimate for the required initial density is overly 
cautious: more accurate estimates which include data on, 
e.g. the emission direction of radiation (our simulations 
show for example that hard photons are typically emit- 
ted away from, rather than toward, the focus '16]) could 
significantly reduce the stated requirements. 

Pair production:- The quantities relevant for pair 
production are not the electric or magnetic field strengths 



themselves, but rather the local Lorentz invariants |4j 
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and V 
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In an e-dipole pulse we have P = E • B = 0, i.e. the 
electric and magnetic fields are orthogonal in all space. 
Near the focus of the pulse, the invariant S is positive 
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A = 1 /im 


A = 0.8 /xm 


A = 0.4 /im 


1660 PW 


1 (5.5 kJ) 


10^ (4.4 kJ) 


10^° (2 k,I) 


1120 PW 


10"'' (3.7 kJ) 


1 (3 kJ) 


10* (1.5 kJ) 


320 PW 


10"^^ (1.1 k,I) 


10""(0.85 kJ) 


1 (0.43 kJ) 



TABLE II. The number of pairs produced (order of magni- 
tude only) as a function of wavelength and power. Bracketed 
figures give the incoming energy in one cycle. Parameters 
are chosen so that the diagonal values exhibit the threshold 
(average) power required to create a single pair in one cycle. 



(the electric field 'dominates') and the pulse is capable 
of pair production O HI [3S] . In general, the key to pro- 
ducing a significant number of pairs is maximising the 
relevant invariant. This can be done by maximising the 
electric field and eliminating the magnetic field in the fo- 
cus, which in practice is achieved (as above) by the use 
of counter-propagating pulses [12]. Define a second 

focusing parameter, f^, via e = ftE^x- "^^^^ relates in- 
coming power to peak invariant e. For two NF beams one 
finds /<: = SttA ~ 2.51 for A = 0.1 (a single NF beam has 
/c — IGttA^ due to a nonzero focal magnetic field), while 
the e-dipole pulse has — Je — ^ — 14.5. Again, the 
dipole pulse gives the most efficient 'conversion' of in- 
put power into invariant e, and so we expect the e-dipole 
pulse to be capable of producing significantly more pairs 
than is possible in other pulses. 

To calculate the number of produced pairs, we use the 
locally-constant-field estimate [TU], based on [23]. For a 
field in which 7^ = we have 



pairs 



47r3A 



(9) 



where e = y^S + \S\/Es, and Ac — h/mc is the (re- 
duced) Compton wavelength. Since the main contribu- 
tion to A^pairs comes from the focal region (which in a 
quasi-Gaussian is localised both in space and time), and 
in order to compare with other pulse models (which are 
typically monochromatic in the far-field), we temporarily 
drop the envelope in ([4]), setting D = oo, and quote only 
the contribution to A'^pairs from a single cycle. 

The number of pairs created is shown in Table [IT] as 
a function of the driving field's wavelength and (aver- 
age) input power. Our results demonstrate that pair 
production can be observed by using several kJ of total 
energy, at optical wavelength: for a driving pulse with 
A ~ 0.8/im, the threshold power required for pair cre- 
ation is around P = 1120 PW, and the corresponding 
peak electric field is Eq — O.OSEs- We see therefore that 
large numbers of pairs can be produced well below the 
Schwinger limit, compare piTT]. 

All predictions of the pair yield based on ([9| should 
of course be interpreted as a measure of how easy it is 
to produce a pair, given an initial power and frequency, 
rather than a number of pairs. The reason for this is that 



once a pair is produced, other processes can occur and 
either increase or decrease the number of pairs able to 
escape the pulse and be detected. In the former category 
is cascade formation [131 US] , in which the pairs are accel- 
erated and emit hard photons which create further pairs, 
and so on [TTl [T2] . In the latter category are processes 
such as pair annihilation, though the annihilation cross 
section in intense fields is typically much smaller than 
that of cascade-generating processes [221 HZ] . It is there- 
fore important to address how one might best detect pair 
production events. We turn to this now, by analysing the 
behaviour of the created pairs. 

Post-creation behaviour:- Fig.[2]shows typical trajec- 
tories, and X values, for positrons born at rest near the 
focus of a 13 = 30 fs e-dipole pulse. Due to the intri- 
cate field structure, the trajectories depend sensitively on 
where the particles are created. The majority of parti- 
cles, with a typical trajectory shown in red in Fig. [2] have 
large x-values, implying that the depicted motion will re- 
ceive significant quantum corrections. We can see though 
that if the particles leave the focus, their x drops quickly 
to below one and classical predictions become accurate. 
Our simulations show that radiation reaction is responsi- 
ble for recirculating particles back into the focus, and has 
an added effect of keeping the particles 'in phase' with the 
field intensity gradient, which is the cause of their large 
X values. Recirculation returns the particles to the quan- 
tum regime, increasing the possibility of additional QED 
processes occurring, as predicted in ^3j. Hence, even 
at intensities well below the Sauter-Schwinger limit, the 
dipole pulse becomes a testing ground not for classical 
electrodynamics, but for QED. 

The behaviour of particles born very close to the fo- 
cus is different. Those particles simply exit the pulse, 
without oscillation, moving almost parallel to the z-axis. 
Further, our simulations show that these particles are 
'out of phase' with the field, so that while they have high 
7 factors, they also have x < !> and so can be analysed 
classically. These particles will emit radiation in a small 
cone about the z-axis. We therefore suggest the detec- 
tion of particles exiting in the ^-direction as a possible 
candidate for direct measurement of the produced pairs. 

Conclusions:- We have considered nonperturbative 
pair production in the fields of an e-dipole pulse. For 
sufficient levels of vacuum, our estimates (which may 
be considered cautious, given the inevitable improve- 
ments in vacuum technology which will be commensurate 
with laser development) show that stray particles are not 
driven into the most intense part of the field, thus reduc- 
ing the chances of cascade initiating and beam depletion. 
This allows higher focal intensities to be reached, and 
allows us to consider the possibility of experimentally 
measuring nonperturbative pair creation, for which a to- 
tal power of 1000 PW is estimated to be necessary. This 
corresponds to field-strengths well below the Schwinger 
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FIG. 2. Trajectories [left] and x parameter [right] for par- 
ticles created in a 30 fs, 1100 PW pulse. Initial positions: 
X = y — z = lO^^A (red), x — y — z = 10~^A (green), 
x = y^z^ lO^^A (blue). 
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